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We investigate multipartite states in the Fenna-Matthews-Olson (FMO) pigment-protein complex 
of the green sulfur bacteria using a Lorentzian spectral density of the phonon reservoir fitted with 
typical parameter estimates of the species, P. aestuarii. The evolution of the entanglement mea- 
sure of the excitonic W qubit states is evaluated in the picosecond time range, showing increased 
revivals in the non-Markovian regime. Similar trends are observed in the evolution dynamics of the 
Meyer- Wallach measure of the JV-exciton multipartite state, with results showing that multipartite 
entanglement can last from 0.5 to 1 ps, between the Bchls of the FMO complex. The teleportation 
and quantum information splitting fidelities associated with the GHZ and Wa resource states of 
the excitonic qubit channels of the FMO complex show that revivals in fidelities increase with the 
J_( ' degree of non-Markovian strength of the decoherent environment. Results indicate that quantum 

information processing tasks involving teleportation followed by the decodification process involving 
Instates of the FMO complex, may play a critical role during coherent oscillations at physiological 
temperatures. 
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I. INTRODUCTION 



During photosynthesis, the vital processes of light absorption, charge separation and efficient energy transfer to 
a reaction center (RC) occurs and culminates in the conversion of electronic energy to useful chemical energy 0- 
6|. These processes are performed by light harvesting (LH) complexes in photosynthetic systems. The well known 
purple photosynthetic bacteria constitutes two different types of LH system units with ring-like structures, the core 
LH1 complex which surrounds each RC, and the peripheral LH2 complexes that transfer energy to the LH1 complex 
043 ■ While the LH1 complex has only one absorption band at around 875 nm, the LH2 complex may have two 
absorption bands (800 and 850 nm) in the case of Rhodobacter (Rb.) sphaeroides. In the Fenna-Matthews-Olson 
(FMO) pigment-protein complex of the green sulfur bacteria, sunlight is transported in the form of electronic energy 
by bacteriochlorophyll (BChl) molecules to the RC complex |10l - fl5| . Light energy is harnessed with a near efficiencies 
of 95% [lj| or more in natural photosynthetic systems, such efficiencies have yet to be realized in artificial light 
harvesting complexes. The observation of long lived coherences, lasting several picoseconds, in photosynthetic antenna 
complexes fl3l Il7l - [25| has been linked to quantum entanglement and elements of quantum communication protocols 
during exciton (correlated electron-hole pair) propagation. 

In general, the exciton propagation times between the photosynthetic pigment molecules are about 100 fs. At room 
• • ■ temperatures, biological systems undergo large decoherences due to interactions with the surrounding environment, 
and excitonic coherences are not expected to persist at the comparatively longer times noted in the experimental 
results. Partly for this reason, there is theoretical interest in the accurate modeling of the dynamics of the exciton 
to examine the underlying basis for efficient energy transfer in photosynthetic structures. However, the energy 
scales involved during the energy derealization and decoherence due to lattice vibrations are almost equivalent 
(10 — 100 cm -1 ), which makes it a challenging task to model photosynthetic systems. The Markovian approximation 
in which an infinitely short correlation reservoir time is assumed, provides convenience in computational anal ysis, but 
yields reliable results only under certain limiting conditions. Consequently, the exciton hopping model 0, H|| [27] has 
been replaced by sophisticated quantum mechanical models of exciton propagation incorporating quantum coherence 
features in recent works [28l - l34j . 

While the coupling of the excitonic states to a dissipative environment is expected to lead to fast decoherence, 
several works have shown that the environment noise can enhance the propagation of energy in light harvesting 
systems [28|, [35l - l40| . Studies [39l - l42T | have shown that attributes such as the non-Markovian interactions and quantum 
correlations present in the phonon reservoir help improve the performance of the light harvesting complex. In Ref. [39] , 
Caruso et. al. showed using analysis of the entanglement properties in the Fenna-Matthew-Olson complex, that the 
delicate interplay of quantum mechanical features and the environmental noise plays a critical role in bringing about 
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an optimal system performance. The authors examined entanglement in further detail in the FMO complex, under 
different natural and artificial conditions and Markovian and non-Markovian dynamics [40| . and showed that near 
unit energy transfer occurs at the optimal system operating point at which the entanglement measure reaches only 
intermediate values. In this regard, a thorough examination of the structure and interactions of quantum states 
in the abstract vector space known as the Hilbert space will contribute to greater insight of optimal processes in 
photosynthetic systems. 

A significant feature of quantum states is their linear structure in Hilbert space. Unusual consequences arise when 
each quantum state can exist as a sum of two distinctly different states in the tensor product space of the two 
separate Hilbert subspaces. This unique property underpins many puzzling properties of quantum states, including 
quantum entanglement which is a valuable resource for the implementation of quantum computation and quantum 
communication protocols (43l - l45j . which include quantum teleportation [46T - I481 ] . dense coding [4{|, quantum cryptog- 
raphy [50l.l5l| and remote state preparation [52J. Moreover, entanglement effects underpins the emergence of classical 
macroscopic features in a system undergoing quantum mechanical interactions within itself and with its surrounding 
environment. The degree of bipartite entanglement of a quantum system which is composed of two subsystems is 
maximum (minimum) when the constituent subsystems are in completely mixed states (pure states). The well known 
concurrence measure |53j is zero for separable states, and becomes one when maximal ignorance of each qubit state 
is achieved in the case of maximally entangled states. A non classical correlation measure known as the quan tum 
discord |54l - l56j was shown to be non-zero even in non entangled states, and more recently, its q-entropies |57j | was 
shown to incorporate aspects of non-classicality not visible within the conventional definition of the quantum discord. 

The exciton manifest as delocalized excitations over the real crystal space [H], l58l - [62| and thus provides the ideal 
example of an extended entangled system in molecular systems [63[. To this end, the exciton can be considered to 
play a central role in the entanglement properties of photosynthetic systems by being in a state of existence at several 
lattice sites and traversing multiple paths simultaneously. The degree of exciton delocalization is however limited 
by the fragility of the coupled molecular chromophores interacting with a bath of environmental phonons and static 
disorder. While the underlying feature that an exciton with a wavevector K is delocalized in real space, is linked 
with the robustness in quantum coherence during solar energy related biochemical reactions, the exact details of such 
investigations has yet to be firmly established in current research on exciton propagation in photosynthetic systems . 
The biggest obstacle to understanding these quantum features is that till now, it is not immediately clear as to how 
entangled excitonic states act to maintain quantum coherence in photosynthetic processes. 

It is important to note the distinction between entanglement measures and those associated with coherences in 
entangled states which evolve under the influence of environmental variables [64-68]. A global system constituting 
several subsystem possess entanglement which decays in ways which are vastly different from the coherence of each 
subsystem, when immersed in a decoherent environment. This stems from the fact that the global characteristics 
of entangled systems and their evolution under local decoherence differ from those of its constituent subsystems. 
Entanglement may vanish at a certain finite time prior to the decay of the coherences [64T4671 169| , and is a topic of 
interest as the use of entanglement as a resource is meaningless if it vanishes rapidly. The abrupt disappearance of 
entanglement at a finite time has been termed entanglement sudden death by Yu et. al. 64]. In general, this finite 
time is expected to decrease for multipartite states or iV-qubit states (N > 2) where several particles are entangled. 
The relationship between coherence and entanglement remains subtle as decoherence in open quantum systems takes 
an infinite amount of time to vanish unlike entanglement, however the degree of robustness of the initial state of the 
quantum system under study is a critical factor that could link entanglement and coherence properties. For instance, 
it is highly likely that robustly entangled states as initial states are more resistent to decoherence processes. 

Unlike bipartite states, investigations related to the robustness of quantum states are challenging for multipartite 
states, where there is still lack of consensus on the categorization and specification of entanglement itself. Recent 
works [7(1 [Zl| have shown variations in the robustness of different types of multipartite states (Dicke, GHZ, W and 
cluster states) under decoherence processes. To this end, there has been scant investigations of the characteristics of 
multipartite states in photosynthetic systems, specifically with respect to the influence of the environmental bath on 
multipartite excitonic states. The study of multipartite states is particularly relevant for light harvesting systems, 
as these states possess a richer source of local and nonlocal correlations due to the multitude of partitions available 
within a group of entangled qubits. Moreover there is relevance in the realistic situation of an entire photosynthetic 
membrane constituting many FMO complexes and thousands of bacteriochorophylls, possibly giving rise to large 
cluster of massively entangled excitonic qubits. In this regard, tt would be worthwhile to seek a detailed understanding 
of how quantum communication protocols (43l . [44| such as quantum teleportation (46l . l47j , are utilized by multipartite 
states to assist energy transfer in photosynthetic processes. The current work is aimed at addressing some of these 
issues. 

In this work, we consider multipartite entanglement in photosynthetic systems and the decay of these states under 
decoherence associated with a phonon bath. A characterization of the different types of multipartite entanglement 
will be examined, including a brief discussion of the robustness of the different types of entanglement against external 
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noise in open quantum system. The influence of specific tasks such as teleportation [46l . |47| . |72| and quantum state 
splitting [73j, [74{ in a noisy photosynthetic system will be examined as well. This work is organized as follows. In Sec. 
ITT1 we discuss various types of multipartite states and the entanglement quantifiers that are used to characterize some 
well known multipartite states. In particular, we discuss in detail the global multipartite state entanglement measure, 
which will be subsequently used to evaluate entanglement quantifiers specific to the FMO complex of P. aestuarii in 
Sec. HID A description of the photophysics and reaction mechanisms in the Fenna-Matthews-Olson (FMO) complex 
from the green sulfur bacteria of the two species, Prosthecochloris (P.) aestuarii and Chlorobium (C.) tepidum is 
provided in Sec. Mil In this Section, we also include a detailed evaluation of the exciton-phonon dynamics based 
on the time-convolutionless (TCL) projection operator technique and a Lorentzian spectral density of the phonon 
reservoir associated with P. aestuarii. Numerical results of the entanglement measure of the excitonic W states of the 
FMO complex species P. aestuarii are presented in Sec. IIII Al for both the Markovian and non-Markovian regimes. 
A similar analysis is performed with respect to the Meyer- Wallach measure of excitonic qubits in the FMO complex 
in Sec. IIVI In Sec. |Vl) . the teleportation and quantum information splitting fidelities of excitonic qubits associated 
with GHZ and Wa resource states, are calculated for the FMO complex species P. aestuarii and the conclusion is 
provided in Sec. (|VI|) . 



II. MULTIPARTITE STATES AND ENTANGLEMENT MEASURES 



A qubit state W associated with N > 2 subsystems can be written in the multi-qubit state or multipartite state as 

N 

I*) = Cn (JO)^™- 1 ' ® |1) ® 10)®^-")) (1) 
n=l 

with coefficients c„. Multipartite systems generally involve the class of permutation-symmetric states which are 
invariant when any pair of particles are swapped, and include the Greenberger-Horne-Zeilinger (GHZ) states [75|-[z3] 
and Dicke states [78| . For a system of N qubits, the Dicke states are obtained using the permutations of basis states 
with N — K qubits being in the |0) state and K in the |1) state using the notation \Sn.k) 

l^ = U) E |0)®|0)®-HQ) |i)®|i)---|i) , (2) 

where < K < N. A well known example of a symmetric Dicke state with just one excitation is the N-qubit W state 

M 

\W N ) = -L (|100...0) + |01...0) + ... + |0...01)) , (3) 



fN 

The W states are inequivalent to the GHZ states [67], [lE Il3 which are obtained as 

\GHZ N ) = a|0)®|0)®..-|0)+^|l>O|l)®."|l) (4) 

S v ' S * ' 

N N 

= a\0f N +P\lf N 

where |a| 2 + \f3\ 2 = I. States derived from each other via local operations and classical communication (LOCC) 
belong to the same group of equivalent resources for quantum information tasks such as quantum teleportation (43[ . 
Consequently, there are differences in robustness against decoherence between entangled W and GHZ states. The W 
state is known to be highly robust in its genuine multipartite entanglement with respect to loss of a single excitation 
[79j . However it is still not conclusive whether the W state is more robust against decoherence than the GHZ states, 
mainly due to the incompatibility in comparison of the dynamics of the two types of states during decoherence. It 
has been shown that in the case of GHZ states undergoing decoherence (68|, the entanglement decay occurs faster 
when the number of initially entangled particles increases. 

Bipartite entanglement are generally quantified by measures such as the von Neumann entropy [13, IsH, negativity 
[H,[I2]' concurrence [53|] and quantum discord [54H56| . These measures cannot be extended to the case of multipartite 
system without introducing an element of subjectivity associated with the partitioning of the global system into smaller 
systems. As a consequence different entanglement measures are expected to display variations in their decay under 
decoherence. Currently, there exists several entanglement measures, such as the the Meyer and Wallach measure 
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and its extension based on the normalized negativity which are used to quantify multipartite entanglement. In a 
recent work, Chaves et. al. [lij identified entanglement measures of states undergoing decoherence with tasks such 
as teleportation H3| and the splitting [73, UM 01 quantum information. Thus in this elegant approach, a task 
associated with quantum information processing was associated with the entanglement quantifiers. 

A commonly adopted approach that is used to determine the multipartite entanglement measure of a 7V-qubit state 
first involves partitioning it 2 N_1 — 1 times. Subsequently, a bipartite entanglement measure such as the von Neumann 
entropy [80l [8lj or negativity (45j | can be evaluated for each bi-partitioned density matrix. The average of the sum 
of the entropies or negativities over all possible partitions may be used to characterize the global entanglement. In 
connection with this, we adopt the global multipartite state entanglement measure used in Ref.[7l| 

, [N/2] 

E = — — V E^ m \ (5) 
\N 72 ^ w 

N% 

E(m) = ( fi ) 

where E(i) denotes the entanglement associated with a single bipartition of the 7V-qubit system and E^ denotes 
the average entanglement between subsets of m qubits and the remaining N — m qubits. The average is performed 
over the N^ l) nonequivalent bipartitions, iV b ™ = ( N n ) if n ^ N/2 and N% /2 = |(^ 2 ) if n = N/2. The total number 

of bipartitions is equal to N cuts — Y^iHi* 1 ^li = 2 JV_1 — 1. Due to the occurrence of mixed states, the negativity is 
employed as a suitable measure of bipartite entanglement. The normalized negativity is defined as [§5[ 

i 

where cti is the negative eigenvalue of the partial transpose matrix for m and the remaining N — m bipartition. The 
measures in Eqs. © and (J7J will be employed in the next section to evaluate the entanglement dynamics of the 
excitonic qubits in the FMO complex. 



III. EXCITONIC QUBITS IN THE FMO COMPLEX 

In this Section, we examine the exciton entanglement dynamics of the Fenna-Matthews-Olson (FMO) complex from 
the green sulfur bacteria of the two species, Prosthecochloris (P.) aestuarii and Chlorobium (C.) tepidum based on 
semiempirical values of the reservoir characteristics at cryogenic temperatures [Hj]. Fig. [T] shows a simplified light 
reaction scheme in which photons are absorbed by a network of antenna pigments linked to the FMO complex, and 
from which the excitation energy is transported to the reaction center (RC) where photochemical reactions convert 
the excitation energy into chemical energy. The FMO complex trimer, which constitutes three identical monomer 
subunits with each unit constituting seven bacteriochlorophyll (BChl)a molecules, acts as an efficient channel of energy 
propagation from the light harvesting molecular complex to the reaction center. 

In the experimental work |86| . FMO complexes from Prosthecochloris (P.) aestuarii and Chlorobium (C.) tepidum 
were isolated and dissolved in suitable solvents before time-resolved transient absorption measurements were per- 
formed. The site energies and coupling between each BChla molecules of Prosthecochloris (P.) aestuarii and Chloro- 
bium (C.) tepidum were characterized Ref. 86 88] using an effective exciton model |87j . which accounts well for the 
optical steady state spectra and the dynamics of the FMO complexes. For both species of FMO complexes, the third 
BChl (numbered 3 in Fig. [1) was found to have the lowest site energy (86|. However there are subtle differences 
between the two species as the absorption spectrum of P. aestuarii has an absorption spectrum maximum height 
at 815 nm, while that of the C. tepidum has its maximum at 806 nm. In P. aestuarii, BChls 1 and 4 contribute 
to the 815 nm absorption spectrum band whereas the contribution to the same band arises mainly from BChl 4 in 
C. tepidum. Consequently, we expect differences in the coherence features that can arise when these two species are 
coupled to a phonon bath with similar environmental parameters. 

In Table U the Bchl site energies for the FMO complex of P. aestuarii (cm ) given in Refs. 
provided. The difference in site energies (provided within { }) is calculated via subtraction from the site energy of 
the third BChl. It is important to distinguish the propagating exciton from the presence of excitation at a specific 
Bchl site, as it is well known that an exciton with a wavevector K is delocalized in real space. The qubit state 
has a probabilistic occupation at all the BChl sites, with one or two BChl sites contributing predominantly to the 
qubit state. We use an effective exciton Hamiltonian defined in the basis of the Q y one-exciton states in which the 
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FIG. 1: Simplified scheme of the light reaction in which photons are absorbed by a network of antenna pigments linked to the 
FMO complex, and from which the excitation energy is transported to the reaction center (RC). Photochemical reactions take 
place at the RC to convert the excitation energy into chemical energy. 



TABLE I: The Bchl site energies (cm -1 ) used in Refs. [Tj (provided as true values of the site energies in last column of Table 
3), [lH and [8(| for the FMO complex of P. aestuarn. The difference in site energies (provided within { }) is calculated by 
considering the site energy of the third BChl as a reference point. 



BChl Exciton energy (cm 1 ) [86J Exciton energy (cm 1 ) [14] Exciton energy (cm 1 ) \15} 

3 12,112 {0 } 12,210 {0} 12,175 {0} 

1 12,266 {154} 12,450 {240} 12,315 {140} 

4 12,293 {181} 12,320 {110} 12,405 {230} 

6 12,396 {284} 12,540 {330} 12,430 {255} 

7 12,457 {345} 12,470 {260} 12,450 {275} 

2 12,496 {384} 12,520 {310} 12,500 {325} 

5 12,634 {522} 12,550 {340} 12,625 {450} 



Q y bandwidth is associated with the lowest excited state of the BChl molecule. When a photon is absorbed by the 
BChl molecule, the newly created exciton propagates rapidly to the adjacent units while interacting with a bath 
of phonons. The exciton can be considered as delocalized over all the seven BChls so that the excitonic states are 
obtained as linear combinations of the excited state wave functions of the individual BChl. In this regard, a simple 
way to represent excitonic states as qubit states is by associating a qubit with the presence (or absence) of an exciton 
at a specific energy level. Each energy level can be modeled as a two-level system, with |0) e (|l) e ) denoting the 

absence (presence) of an excitation at the specified energy level. 

In order to compute the contributions from different pigments to the seven exciton states, we model the subunit of 
the FMO complex via the Hamiltonian, H ex in the site basis, with the coupling energy terms given in units of cm -1 
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The off-diagonal terms in H ex are obtained using the intersite excitonic couplings between various BChls as given 
in Ref.fTij for the FMO complex of P. aestuarii. The various coupling energies were evaluated [Tij in a dielectric 
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TABLE II: Delocalized exciton qubit states provided as linear combinations of probability occupation amplitudes associated 
with the seven BChl sites. The excitonic qubit states are labeled according to increasing exciton energies for the FMO complex 
of P. aestuarii. 
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environment by solving a Poisson equation for each BChl by a finite difference method. The diagonal terms of H ex in 
Eq. ((8]) were obtained using the energy differences for the respective BChl site using the site energies, provided within 
{ }, in Table |H The exciton energies and eigenstates of H ex in Eq. (j8]) can be evaluated using standard diagonalization 
techniques. The seven excitonic qubits (labelled in order of increasing energy) are given in Table [TT] It can be seen 
that each exciton state is generally associated with 1-3 Bchl sites, and the lowest energy exciton qubit state |l) e 

with energy —24 cm -1 is localized mostly on BChls 3 and 4. Most importantly, Table [TTI shows that each exciton 
qubit state may have simultaneous occupation of two or more BChl sites, a feature that is only feasible in quantum 
systems. 

In general, the dimer model consisting of pair of spin-boson [89j j can be used to simplify analysis of the entanglement 
dynamics of the excitonic qubits. The spin-boson system is a well known quantum dissipative system which constitutes 
two energy levels that couple to an infinite system of harmonic oscillators. The case of the multipartite entanglement 
can be examined by extending the global system to several spin chains. We model each two-level excitonic qubit 
interacting with phonons via the Hamiltonian, 

H = woo-+o-- X] hu <i fo q fo q + X! K ( a - & q + a + & q) ' ( 9 ) 
q q 

where a + = |1)(0| and er_ = |0)(1|, are the respective Pauli raising and lowering operators of the exciton with 
transition frequency, ujq. 

In order to utilize Eqs. © and ([7]) for further analysis of the excitonic qubits, and to keep the problem tractable, 
we assume that each qubit is coupled to its own reservoir of phonons. The phonon reservoir is given by the second 
term on the right hand side of Eq. ([9]) where 6 q and & q are the respective phonon creation and annihilation operators 
with wave vector q. The last term of Eq. ^ denotes the qubit-oscillator interaction Hamiltonian which we assume 
to be linear in terms of the phonon operators. A q is the coupling between the qubit and the environment and is 
characterized by the spectral density function, J(uS) — A^<5(w — w q ). 

For the spectral density associated with the coupling of the BChl exciton to the phonon bath environment, we 
consider a Lorentzian spectral density of the reservoir fitted with estimates derived from site-selective fluorescence 
measurements on the 825 nm band of the FMO complex of P. aestuarii [86j j 

m = ( 2 j (io) 

whose central peak is detuned from the exciton transition frequency too by an amount 8 which is treated as a parameter 
for subsequent analysis. It important to distinguish between the molecular transition frequency associated with the 
Q y band of the Bchl-a molecular system (about 12100cm- 1 , see TableU) and the actual excitonic energies which are 
of the order of the typical binding energies of 0.1-0.2 eV in organic molecules [59(. We consider that the phonon modes 
have maximum energies just enough to convert the exciton into free electron- hole pairs, however in general the greater 
involvement of a lower range of phonon energies of the order 0.01 eV (or 80 cm -1 ) is assumed. As a standard, we use 
the full width at half- maximum Aw=80 cm -1 as used in Ref.[86[, noting that the width is related to the reservoir 
correlation time tb via Au = -2- . However in order to examine the influence of the reservoir correlation time tb on 
population differences, we parameterize Aw around the region where phonon energies lie close to 80 cm -1 . 
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The parameter 70 is associated with the relaxation time scale tr of the exciton via the relation tr — Jq 1 , and 
7o=lllps~ 1 was used in Ref. 86]. 70, is expected to increase with temperature, hence for subsequent calculations, 
this term will be treated as a parameter. While the conditions under which the experimental work in Ref. (8(| was 
carried out precludes significant temperature-dependent dephasing processes, we consider higher values of 70 (1000 to 
2000 cm^ 1 ) to incorporate non-Markovian effects. Hence we assume that the spectral density profile in Eq. ([9]) holds 
valid at higher temperatures (> 50 K). 

Due to the exciton-phonon interactions, the excitonic qubit decays to oscillator states in the reservoir, making a 
transition from its excited state |1) to ground state |0) . We assume an initial state of the qubit with its corresponding 
reservoir in the vacuum state of the form [9C| 

JV' 

|0i> = ]l> e (S> JJ I0 fc > r = |l> e cs> ]0> r , (If) 

k=l 

where |0) implies that all N' wavevector modes of the reservoir are unoccupied in the initial state. The subscripts 
e and r refer to the excitonic qubit and the corresponding reservoir respectively. \<pi) undergoes subsequent decay of 
the following form 

|&>— >u(t) \l)JO) I + v(t) |0)Jf) r , (12) 
In order to keep the problem tractable we consider that |1) denotes a collective state of the reservoir as follows 

i 1 >' = iE A wWiw), (13) 

^ ' n 

where {n} denotes an occupation scheme in which there are oscillators with wavevector k = i in the reservoir and 
we define the state \{n}) as 90] 

|{n}) = |no 1 n ll n 2 ...n i ..niV') ) (14) 

In the collective reservoir state, the phonon oscillators can be present at all allowed modes, including simultaneous 
excitation of several phonon states. At non-zero temperatures, the reservoir state is a Boltzmann weighted average 
over all possible permutations of the occupation scheme {n}. 

Following the time-convolutionless (TCL) projection operator technique described in Chapter 10 of Ref.[9l[, we 
obtain an expression for the coefficient u(t) (see Eq. (fTTj)]). 

t 

u(t) = — J J cLjJ(uj) exp[i(uJo — uj)(t — ti)] u(ti)dti, (15) 


which can be solved by means of a Laplace transformation. The time-convolutionless projection operator technique is 
independent of the choice of the spectral density of the environment and accounts for strong non-Markovian dynamics 
[9l| . that may arise in the photosynthetic molecular environment. 
Eq. (|15p is further simplified using the form of spectral density in Eq. (fT0|) 

t 

u(t) =-[ i 7 o(A^/2 - i5 ) e -{^/2-i6)\t-t x \ u{ t l)dtlj (16) 


The Laplace transformation of Eq. (fl"5)) with initial condition u(0) = 1, yields the solution [9l[ 



u(t) = e -(A«/2-tf)t/2) 



£A Aw/2 - iS . (it 
cosh — M : sinn ' 



2 J £ V 2 



(17) 



where ^ = yj (Aw/2 — id) 2 — 70 Aw. Fig.[2ji,b,c illustrates the evolution of the coherence properties as reflected in the 
exciton population difference, AP between the excited state |1) and ground state |0) , evaluated using |u(t)| 2 (Eq. 
(fT7| . The change of AP with time t (ps), Aw/2 and 70 shows that for smaller Aw/2 ~20 cm -1 (or larger reservoir 
correlation time) and large 70 (or small exciton relaxation times), there is increased time period (up to 1 ps) over 
which the population difference, AP which is a signature of coherence, is maintained for a select range of parameters. 
The time duration of effective coherence (about 1 ps) obtained here appear to compare well with times of population 
transfer between various exciton levels (about 1 ps for levels 1 and 2) noted for the FMO complex of P. aestuarii 
[86] ]. These results suggest that the non-Markovian regime is best suited for occurrence of large coherence times. The 
analytical form for u(t) in Eq. (| 17|) will be useful for further analysis of the entanglement measure for the W states 
of the FMO complex which we describe next. 
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(a) AP 




FIG. 2: (a) The population difference, AP between the excited state |l) e and ground state |0) e , evaluated using u{i) 2 (Eq. 
(|17|) . as a function of time t (ps) and Aoj/2 (cm -1 ) at 70 =1000 cm" 1 and detuning parameter, S =0. (b) AP as a function 
of time t (ps) and 70 (cm" 1 ) at Auj/2 =40 cm" 1 . (c) AP as a function of time t (ps) and 70 (cm ) at Aw/2 =20 cm" 1 . 



A. Entanglement dynamics of excitonic W states of the FMO complex 



We consider a system of four non-interacting excitonic qubits coupled to uncorrelated reservoirs and which are 
initially prepared in the W state as follows 

\M = (|0001) ei e 2 e 3e4 +|0010) ei e 2 e 3e4 + |0100) eie2e3e4 (18) 

+ |1000) ei e 2 e 3 e 4 )|0000) rir2r . 3r4 /2, 

The system of four qubit state is a reasonable compromise between a high number of entangled states and resource 
available for numerical computations. The evolution of the total system can be obtained using Eq. (1121) . so that a 
typical form of the following term evolves as 

|0)e 2 I0)e 3 I0)e 4 |0> r , (O^ 10)^ , 

(19) 

where u\(t) and v\(t) (see Eqs. (|12[) . (JT7J)) are interaction parameters associated with the first excitonic qubit which 
we label via the subscript 1. Continuing with the same manner with the remaining three excitonic qubits, we obtain 

\fo) = [|000) eie2e3 |000> rir2r3 ( U4 (t)|l>e 4 |0) r4 + W 4(t)|0)e 4 |l>r 4 ) 

+ l° 00 >e ie2e4 I000) rirar4 («a(*)|l>« 8 |0), 3 ) + «3(t)|0)e, |l)r 3 ) 
+ l° 00 )e ie 3e 4 I000)r 1 r s r ll M«)|l)e 2 |0) ra ) + U 2 |0)e 2 |l) r2 ) 

+ l° 00 >e,e,e 4 |000> («l (*) 1 1> ^ |0> n + V X (t) |0) e± | 1) n )] /2. (20) 

The reduced density operator of the excitonic qubit subsystem p e (t) = Tr r (\<f>t) (4>t\) and the reduced density operator 
of reservoir subsystem p r (t) = Tr e (|</) t ) (4>t\) can be obtained, provided details of each individual system-reservoir 
parameters, Ui,vi (i=l-4) are known. The detailed derivation of the density operator matrices for such generalized 
exciton-phonon reservoir systems is complex and for ease in numerical evaluations, we consider the simple model in 
which u\ = U2 = u-3 = = u. Thus each of the exciton is assumed to be coupled to a independent phonon bath, with 
same spectral density. This allows convenient evaluation of the entanglement measure of the excitonic qubit (denoted 
by E e ) and phonon reservoir E r systems, using Eqs.© and ([7]), and at the same time is not a significant compromise 
on the accuracy of the evaluated results. Three parameters which appear in the spectral density relation in Eq. (|10[) 
are expected to play critical role in the analysis of the entanglement dynamics, (a) time (t), (b) 70 which is associated 
with the exciton relaxation time and (c) S which is the amount by which the central peak is detuned from the exciton 
transition frequency. 

We consider both the Markovian (Aw ^70) and non-Markovian regimes (Aw <C 70) as well as the non- resonant 
case for which 6^0. Fig. [5Ji,b illustrates the evolution of the entanglement measure of the excitonic qubit subsystem 
(reservoir subsystem), E e (E r ) as function of time t (in picoseconds) and 70 (in units of cm" 1 ). The entanglement 
measures are evaluated for the W state (Eq. (fT5|) using Eqs.©, ([7]) and (IT7|) . and setting the detuning parameter 
at 8 =0. It can be seen that as 70 increases, the sub-systems move into the non-Markovian regimes, and increased 
oscillations in revivals of entanglement, associated with comparatively long phonon reservoir correlation times, can 



|l)e, |0)e. I°)e„ |0)e. |0) r . |0> r . |0) r . |0>, 



lo>,, +«i(*)|o> e . |i>, 
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be seen in both sub-systems. The non-Markovian features are expected due to inflow of quantum information from 
the reservoir subsystem back into the excitonic qubit subsystem. Overall there is net flow of entanglement from the 
excitonic subsystems to the phonon reservoir systems after 0.5 ps. Fig. |3j; also illustrates that for smaller Aw/2 ^20 
cm -1 where there is higher reservoir correlation time, there is a pronounced revival in the entanglement measure E e . 

Figs. |4ji, b also show that both the entanglement measures, E e and E r decrease when the detuning 6 increases, 
which has been attributed to the effective decrease in coupling between the two interacting subsystems [92j . The time 
duration for which significant entanglement exists (~ 0.1-0.3 ps) obtained here is less than the decoherence times 
noted for the FMO complex of P. aestuarii [861 ] . The evaluation of the entanglement measures, E a and E r based on 
GHZ states are expected to yield similar qualitative variations with respect to change in time t, 70 and d as obtained 
for the W states. However some quantitative differences due to the different degree of robustness of the different 
multipartite states, can be expected in the entanglement measures, E a and E r . 




FIG. 3: (a) Entanglement measure of the excitonic qubit subsystem, E e as function of time t (ps) and 70 (cm ), evaluated 
for the W state (Eq. (fl8)) using Eqs.([5| and (0. The detuning parameter, 8 =0. Alj=80 cm -1 [8y] in Eq. (Tl0)| . The parameter 
estimates are typical of the FMO complex of P. aestuarii [86| (b) Same as in (a) except that the entanglement measure of 
the system of the collective reservoir states E r is obtained as a function of t (ps) and 70 (cm -1 ) (c) Entanglement measure 
of the excitonic qubit subsystem, E e as function of time t (ps) and Alj/2 (cm -1 ), with the detuning parameter, 8 —0 and 70= 
1000 cm" 1 . 




FIG. 4: (a) Entanglement measure of the excitonic qubit subsystem, E e as function of time t (ps) and 8 (cm ), evaluated for 
the W state (Eg .1(15)1 using Eqs.j5} and ©. 70 is set at 1000 cm -1 and Ao;/2=40 cm -1 . 

(b) Same as in (a) except that the entanglement measure of the system of the collective reservoir states E r is obtained as a 
function of t (ps) and 8 (cm -1 ) 

It is important to note that the results obtained here are based on simplified model systems with underlying assump- 
tions, such as the requirement that all qubits in different Hilbert subspaces experience similar isotropic interactions 
with the phonon reservoir subsystem. The qubits are assumed to be independent of each other, and coupled to 
uncorrelated reservoirs. In the realistic situation, each excitonic qubit may share one or more BChl sites as shown in 
Table HT1 which means that the reservoir subsystems are correlated as well. Due to the large energy difference between 
each qubit system (about 100 cm -1 , Table [TTJ) , the energy profile of phonon modes linked to different qubit are ex- 
pected to be different, which to some extent justifies the assumption that the reservoir system are uncorrelated. The 
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basic models employed here also circumvent the inclusion of more realistic features such the correlated fluctuations 
of the exciton transition energies. It has been shown that environmental changes lead to correlated fluctuations of 
the exciton transition and interaction energies with impact on population transfer rates, decoherence rates and the 
efficiency of photosynthetic complexes [25[ . Even interaction fluctuations of a small magnitude results in an increase 
in efficiency of the FMO complex (25|. Nevertheless, the use of the assumptions mentioned earlier have enabled the 
numerical evaluation of some estimates of the timescales involved in the entanglement dynamics of the system under 
study, and to help understand on a qualitative level some critical parameters responsible for quantum coherence in 
light-harvesting systems. The incorporation of realistic features will lead to slight changes in the magnitude of the 
entanglement measures considered here and in the next Section, however the qualitative features are expected to be 
retained. 



IV. MEYER- WALLACH MEASURE OF EXCITONIC QUBITS IN THE FMO COMPLEX 



In this section, we evaluate the Meyer- Wallach measure for a system of photosynthetic qubits, with each qubit 
interacting with its own reservoir of phonons. This monotone measure was defined by Meyer and Wallach [83j as a 
single scalar measure of pure state entanglement for a system of more several qubits. For the case of N qubits, the 
Meyer- Wallach measure is based on the entanglement of each qubit with the remaining (N — l)-qubits as follows 



1 N 

-Y2(l 
Nf-t v 



(21) 



k=l 



where is the reduced density matrix of the fcth qubit obtained after tracing out all the remaining qubits. The Meyer- 
Wallach measure has inherent deficiencies in that it is unable to distinguish states which are fully inseparable from 
states which are separable into states of subsystems. Nevertheless, it is a useful qualitative measure in the investigation 
of subtle features of entangled systems, and which we will utilize to analyze the dynamics of the photosynthetic system 
under study in this Section. 

We first consider the joint evolution of a pair of non-interacting excitonic qubits coupled to uncorrelated phonon 
reservoirs, and which are initially prepared in the following state 



bo) = [a(|00) ei e 3 +6|ll) ei e 2 ] |00) P1 , 



(22) 



The state above evolves as a four-qubit multipartite state influenced by the real coefficients a, b. There are obvious 
differences between the entangled states in Eq. ([22|) and (|18p. The density matrix of the excitonic subsystem derived 
from the global state in Eq.(f2"2"j) and expressed within the basis (|0 0) , |0 1) |1 0) |1 1)) appear as 



.(*) 





hit) 
o 








hit) 





(hit) 
o 



\hit) 

For t > 0, the matrix elements of the matrix appearing in Eq, 

h{t) = a 2 + b 2 v 1 (t) 2 v 2 (t) 2 
hit) = b 2 Vl it) 2 U2 (i) 2 , 

hit) 

hit) 



hit) \ 
o 



hit) J 

evolve as 



(23) 



b 2 mit) 2 v 2 it) 2 , 

b 2 Ul it) 2 u 2 (t) 2 , 



hit) = ab ui(t) u 2 (t). 

where Ui,Vi are conversion functions associated with the ith exciton and its phonon reservoir. The problem can be 
extended to the photosynthetic system examined earlier by substituting u with the form obtained in Eg. (1171) . It 
should be noted that the X-state structure in Eq. (|2"3"|) preserve its X-form during evolution. 

For density matrix given in Eq. (123p , an explicit expression of the Meyer- Wallach measure for the iV-exciton multi- 
partite state is obtained as [90l . 193 



Q = 2a 2 b 2 + Ab 2 u 2 v 2 



(24) 



where it is assumed that Ui = u for all the N exciton-reservoir interacting system. Due to this assumption, the 
Meyer- Wallach measure, Q in Eq. (|2"4"|) is independent of N. For the elaborate and more realistic models, in which some 
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exciton-reservoir systems evolve differently from adjacent pairs, Q becomes a complicated function of the differences 
in system-reservoir dynamics in separate Hilbert subspaces. For this reason, we chose a highly simplified scenario of 
the same mode of exciton-reservoir for all N qubits, as was also the case in Sec. MI Al 

Substituting the form of u given in Eq. lfTT)) into Eq.([2"4"|). the evolution dynamics of the Meyer- Wallach Q measure 
is plotted in Fig. [5^i,b. The figure shows expected differences between the Q measure and the entanglement measures 
E e and E r obtained in the preceding Section. Unlike E e , the Q measure incorporates the entanglement of the global 
system which constitutes the reservoir subsystems as well. As a consequence, the magnitude of Q shows slower decay 
with time t at large 70 ~ 800cto _1 as seen in Fig. [5^,. Q is also sensitive to the correlation parameters a and b (see 
Eq. (l22j0 as illustrated in Fig. [5Jd. As in the case of the entanglement measures E e and E r , there is increased revivals 
in entanglement as the degree of non-Markovian strength increases. These revivals can be seen as a signature of 
coherence, and hence the correlation parameter b plays a critical role in the length of time over which such revivals 
can be maintained. In general, the maximum Q = 1 is obtained for a = 0, b — 1 with u\ = U2 = -4?. 




FIG. 5: (a) MeyerWallach measure Q (Eq.([24])) as a function of time t (ps) and 70 (cm -1 ), evaluated for the state in Eq. (|22[) 
using Eq. (I24|l . a=0, 6=1 and Aa>/2=40 cm -1 . (b) MeyerWallach measure Q as a function of time t (ps) and parameter b, 
70=800 cm -1 and Ao;/2=40 cm -1 . (c) Meyer Wallach measure Q as a function of time t (ps) and Aw/2, 7o = 1000 cm -1 , 
a=0, b=l. 



V. TELEPORTATION AND SPLITTING FIDELITIES OF EXCITONIC QUBITS IN THE FMO 

COMPLEX 

In this section, we apply the results of teleportation fidelities obtained by Chaves et. al. [IH to analyze the 
entanglement dynamics of the FMO complex detailed in Sec. Mil Using the concept that multipartite entanglement 
is linked to the teleportation and quantum information splitting protocols [84j . we examine the effect of the phonon 
reservoir on the entanglement properties of the excitonic qubits in the FMO complex. This is also linked to the 
possibility that the entangled qubit states in the photosynthetic complex may act as a resource for the teleportation 
and quantum information splitting protocols. 

In the teleportation process pa - |48j , two parties conventionally referred to as Alice and Bob share an entangled 
bipartite state, with the quality of the teleportation measured by the fidelity / = (ijj\p\ip), where p is the state obtained 
after teleportation of an unknown quantum state \ tp) . For the unknown state of dimension d. the maximum achievable 
fidelity / max for a given bipartition of the state used as the medium of transmission is [84| / max < 2 ~d+i wnere t ne 
negativity M is associated with the given bipartition. The bipartite mixed state leads to imperfect teleportation, and 
thus reduce the fidelity of teleportation to values below those of the pure states which yield a fidelity of § [48| . 

Chaves et. al. considered a viable connection between global entanglement quantifiers in entangled systems and 
the robustness of tasks such as teleportation and the splitting of quantum information [73ll74j|. which is a generalization 
of the teleportation protocol. During the splitting process, an unknown qubit |^o) = cos (9/2) |0) + e 1 ^ sin (0/2) |1) is 
sent to N other parties via a shared entangled state. The no-cloning theorem ensures that only a single copy of |^o) 
can be retrieved. The association between the degree of cooperation of the N parties during the extraction process 
and the global entanglement of the shared GHZ or W-like states under decoherence were examined in Ref. [84( ■ The 
results that is of particular relevance in the context of photosynthetic systems, are the explicit forms of teleportation 
and quantum information splitting fidelities obtained in Ref. |84| for protocols associated with GHZ and Wa states, 
under the amplitude-damping channel. The Wa state is associated with N + 1 qubits and of the form 

(25) 



1 

77? 



|o> \w N ) + |1) |0) 
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We note the similarities between the amplitude-damping channel and the action of the phonon reservoir on the 
photosynthetic qubit as detailed in Sec. IIIII Within this channel, the population of the excited state is reduced by a 
factor 1 — p, where the p is a parameterization of time. For the excitonic qubit interacting with a phonon bath, one 
has p = 1 — \u(t)\ 2 , where u(t) appears in Eq. p7|) . At initial time t = 0, u(t) = 1 and p = 0, and at t— > oo, u(t) = 
and p = 1 . 

The teleportation fidelity, averaged over all the possible input states l^o), was obtained as 



F, 



1 



<•'//:. g 



(l-pf- 1 (2-p) + 2(l-p) N/2 +p N - 1 (l+p) (26) 



for GHZ resource states. In the case of the Wa state associated with N + 1 qubits, = 
)0)\W N ) + \1)\Q)* 



i 

V2 



, the teleportation fidelity was obtained as [8 

1 



F w =~(3-2p + p 2 ) (27) 

Chaves, et. al. [Hf likewise obtained simple form for the fidelities associated with the complete splitting protocol 
involving teleportation followed by the dccodification as 



F = I 

GHZS ^ 



2-p(l-p) + (l-p) N/2 (28) 



F ws = (29) 

The full details of the derivation of the fidelities in Eqs. (gH), (EH © are provided in Ref.[s3|. 

The evolution dynamics of the teleportation fidelities for the systems under study here, can be obtained as shown 
in Fig. [5] and [5^ by substituting the form of u (Eq. flTTl) ') into Eqs. ([2l)|) and (|2"7)l . and using typical parameter estimates 
of the FMO complex of P. aestuarii [§6| . Likewise the fidelities associated with the complete splitting protocol 
(teleportation followed by the decodification) obtained using Eqs. flTTj) . (|2"8|) and (|2T))) are plotted in Fig. [7] and [5}}. 
Both the GHZ and Wa resource states reach the classical fidelity of 2/3 in the Markovian and non-Markovian regimes, 
for the two types of fidelities at t— > oo [84j . The teleportation fidelities associated with the Wa resource states lie 
above the classical threshold fidelity of 2/3, while those of the GHZ resource states lie below the classical threshold 
fidelity. The results obtained here are in line with earlier observations by Chaves et. al. [H| that the Wa states 
display greater robustness than the GHZ states. 

The increased revivals in the fidelities with time, for increasing degree of non-Markovian strength, can also be 
noted for both GHZ and Wa states, with the amplitude of oscillations bounded either from above or below by the 
classical threshold value. These revivals, appear to last up to time periods of typically 0.5 - 0.7 ps, depending on 70, 
which is associated with the exciton relaxation time. The oscillations in the fidelity measure is indicative of strong 
interactions between the excitonic qubit and the phonon reservoir system. Fig.[6K,b,c shows the overall decrease of the 
teleportation fidelities with the number of qubits N in the GHZ state, however there is also an associated increased 
oscillations when N is increased at the non-Markovian regime. Unlike the teleportation fidelities, Fig. [7^,b,c shows 
that the splitting fidelity, F GHZS is more robust to changes in the number of qubits TV in the GHZ state. Comparison 
of the results in Figs. [BJ [7] and E] show that in general, higher fidelities are obtained in the case of the complete splitting 
protocol involving teleportation followed by the decodification. In particular, there is notable increased oscillations 
of the fidelities associated with the complete splitting protocol in the non-Markovian regime when 70 is increased to 
1500 cm" 1 (see Fig. |5}. 

Any connection between coherence and fidelity appears to be subtle, and a rigorous link between the two measures 
have been not demonstrated so far, and is beyond the scope of this work. As the fidelity characterizes the quality 
of information transmission through the appropriate quantum channels (e.g. excitonic qubit channels of the FMO 
complex), it would not be naive to link higher fidelities with more robust coherence phenomenon. Using this reasoning, 
we can attribute coherent oscillations at physiological temperatures to quantum information processing tasks involving 
teleportation followed by the decodification process involving Wa states of the FMO complex. Such tasks in particular, 
may become more important at large 70, which is more likely to occur at higher temperatures. 

It well known that one antenna complex serves many FMO complexes, and multiple FMOs are connected to a given 
reaction center, a structural arrangement which allow photosynthetic organisms to most effectively organize cellular 
resources @, 0, 12H . Even though the inter- monomer coupling strengths are about one order of magnitude smaller 



than intra- monomer couplings 14|, the three weakly coupled subunits in the FMO photosynthetic complex trimer, 
can form a larger cluster of excitonic qubits. There also exists the occurrence of massive entangled excitonic qubits of 
N > 64 across an entire photosynthetic membrane constituting thousands of bacteriochorophylls as determined using 
atomic force microscopy in Ref. |94j |. 
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In summary, there are increased revivals in fidelities in the non-Markovian regime for resource states (GHZ states, 
IT-like states) in contact with a decoherent environment. This can be attributed to the comparatively long phonon 
reservoir correlation times associated with non-Markovian interactions. Unlike the GHZ states, the IT-like states 
are expected to provide greater resilience to environmental related decoherence during photosynthetic processes. The 
contribution to global entanglement due to teleportation followed by the decodification tasks in Wa states of the FMO 
complex appear to play a domineering role during quantum oscillations at physiological temperatures. The possibility 
of the latter processes accounting for experimental results [TH, [2(| which show that coherent oscillations persist up to 
1 ps at physiological temperatures cannot be ruled out. These results have important implications for multipartite 
states present in the FMO complex of light-harvesting systems. 
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0.6 0.8 
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FIG. 6: (a) Teleportation fidelity, F GHZ , as a function of time t (ps) and 70 [10 cm - (blue line), 50 cm -1 (green line), 1000 
cm -1 (red line)] based on the GHZ resource states, using Eq. (I26[) . The detuning parameter is set at 8 —Q and N=4. Aw=80 
cm -1 ^86j in Eq. (flO)) . The classical fidelity of 2/3 is denoted by dotted lines. (b) Same as in (a) except iV=16 (c) Same 
as in (a) except iV=64 
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FIG. 7: (a) Quantum information splitting fidelity, F GHZS , as a function of time t (ps) and 70 [10 cm" 1 (blue line), 50 cm -1 
(green line), 1000 cm" 1 (red line)] based on the GHZ resource states, using Eq. (I28[) . The detuning parameter is set at 8 =0 
and 7V=4. Aoj=80 cm -1 [861 ] in Eq. (|10[) . The classical fidelity of 2/3 is denoted by dotted lines. (b) Same as in (a), except 
iV=16 (c) Same as in (a) except A?=64 



VI. CONCLUSION 



In this work, we examined the exciton entanglement dynamics of the Fenna-Matthews-Olson (FMO) pigment- 
protein complex from the green sulfur bacteria of the species, Prosthecochloris (P.) aestuarii using typical values of 
the reservoir characteristics at cryogenic and physiological temperatures [8a ]. Using the time-convolutionless (TCL) 
projection operator technique and a Lorentzian spectral density of phonon reservoir, the evolution of the entanglement 
measure of the excitonic qubit W states is evaluated, showing increased oscillations in the entanglement in the non- 
Markovian regime. The calculations are repeated for the evolution dynamics of the Meyer- Wallach measure of the 
iV-exciton multipartite state showing that multipartite entanglement can last up to at least 0.5 to 1.0 ps, between the 
Bchls of the FMO complex, in a decoherent environment. 

We have also considered quantum information processing based on teleportation and a complete protocol involving 
for the GHZ and Wa resource states, within the context of photosynthetic systems. The teleportation fidelities 
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FIG. 8: (a) Teleportation fidelity, F w , as a function of time t (ps) and 70 [10 cm -1 (blue line), 50 cm -1 (green line), 1500 
cm -1 (red line)] based on the Wa resource states, using Eq. (|27|l . The detuning parameter is set at 5 =0 and Au;=80 cm -1 
[8t3 | in Eq.H10[). The classical fidelity of 2/3 is denoted by dashed lines. 

(b) Quantum information splitting fidelity, F ws , as a function of time t (ps) and 70 [10 cm -1 (blue line), 50 cm -1 (green line), 
1500 cm -1 (red line)] based on the Wa resource states, using Eq. (|29[). The detuning parameter is set at S —0 and Auj—80 
cm" 1 [H| in Eqs. lfTU)) . 



associated with the GHZ and Wa resource states associated with the excitonic qubits in the FMO complex of P. 
aestuarii show increased revivals in the fidelities as the degree of non-Markovian strength is increased. Unlike the 
GHZ states, the U/-like states appear to provide greater resilience to environmental related decoherence during 
photosynthetic processes. Results indicate that quantum information processing involving teleportation followed by 
the decodification tasks in Wa states of the FMO complex is likely to account for experimental results which show 
persistence of coherent oscillations at physiological temperatures. 

The results obtained thus far highlight the importance role of several parameters (Aw associated with the phonon 
reservoir correlation times, 70 associated with the exciton relaxation times, 5 associated with the detuning parameter) 
involved in the long coherence times of multipartite states in the Fenna-Matthews-Olson (FMO) pigment-protein 
complex of the green sulfur bacteria. In particular the non-Markovian regime appears best suited for occurrence of 
large coherence times for the model system used in this study. We note however that processes in non-Markovian 
systems are not well-understood, and further analysis is need to confirm whether non-Markovianity, associated with 
the backfiow of information from the reservoir to the excitonic qubit system, is indeed instrumental in bringing about 
fast energy transport during photosynthesis. 

Lastly, in large light harvesting antennae systems, with intricate network connectivity, coherence between several 
exciton states involving multipartite states can occur and contribute to the unique ability of these systems to attain 
robustness against decoherence. The analysis of multipartite states carried out in this work can be extended to esti- 
mate the multipartite entanglement measures in photosynthetic systems with alternative forms of spectral functions 
associated with the environmental couplings . The results obtained in this work may provide useful guidelines for 
future experimental work involving the detection of multipartite entanglement in photosynthetic systems. In partic- 
ular, advanced measurement tools involving spectrally resolved, 4-wave mixing measurements [95| may be utilized 
to examine the remarkable resilience of quantum correlations, and the role of multipartite states in photosynthetic 
systems in contact with a decoherent environment. 



VII. ACKNOWLEDGMENTS 



This research was undertaken on the NCI National Facility in Canberra, Australia, which is supported by the 
Australian Commonwealth Government. 



[1] R. S. Knox, in Primary Processes 0} Photosynthesis, edited by J. Barber, pp. 5597 (Elsevier, Amsterdam, 1977). 

[2] H. Zuber and R. Cogdell, in Anoxygenic Photosynthetic Bacteria, edited by R. Blankenship, M. Madigan, and C. Bauer, 

pp. 315348 (Kluwer Academic, Dordrecht, 1995). 
[3] B. Chance and M. Nishimura, Proc. Nat. Acad. Sci. U.S.A., 46, 19 (1960). 



15 



[4] Light-Harvesting Antennas in Photosynthesis, edited by B. R. Green and W. W. Parson (Springer, New York, 2003). 

[5] V. May and O. Kiihn, Charge and Energy Transfer Dynamics in Molecular Systems, 2nd ed. (Wiley- VCH, Berlin, 2004). 

[6] H. van Amerongen, L. Valkunas, and R. van Grondelle, Photo synthetic Excitons (World Scientific, Singapore (2000). 

[7] V. Sundstrom, T. Pullerits and R. van Grondelle, J. Phys. Chem. B 103, 2327 (1999). 

[8] G. R. Fleming and G. D. Scholes GD Nature 431, 256 (2004). 

[9] R. J. Cogdell, A. Gall A and J. Kohler, Q. Rev Biophys 39, 227 (2006). 
[10] T. Brixner, J. Stenger, H. M. Vaswani, M. Cho, R. E. Blankenship, and G. R. Fleming, Nature, 434, 625 (2005). 
[11] R. E. Fenna and B. W. Matthews, Nature, 258, 573 (1975). 

[12] A. Camara-Artigas, R. E. Blankenship, J. P. Allen Photosynth. Res. 75, 49 (2003). 

[13] G. S. Engel, T. R. Calhoun, E. L. Read, T.K. Ahn, T. Mancal, Y. C. Cheng, R. E. Blankenship, and G. R. Fleming, 

Nature, 446, 782 (2007). 
[14] J. Adolphs and T. Renger, Biophys. J. 91, 2778 (2006). 

[15] M. Wendling, M. A. Przyjalgowski, D. Giilen, S. I. E. Vulto, T. J. Aartsma, R. van Grondelle, and H. van Amerongen, 

Photosynth. Res. 71, 99 (2002). 
[16] Y. C. Cheng and G. R. Fleming, Annu. Rev. Phys. Chem., 60, 241 (2009). 

[17] M. Cho, H. M. Vaswani, T. Brixner, J. Stenger, and G. R. Fleming, J. Phys. Chem. B 109, 10542 (2005). 

[18] G. Panitchayangkoon, D. Hayes, K. A. Fransted, J. R. Caram, E. Harel, J. Wen, R. Blankenship, G. S. Engel, Proc. Nat. 

Acad. Sci. 107, 12766 (2010). 
[19] H. Lee, Y.-C. Cheng, and G. R. Fleming, Science 316, 1462 (2007). 

[20] E. Collini, C. Y. Wong, K. E. Wilk, P. M. G. Curmi, P. Brumer, and G. D. Scholes, Nature, 463, 644 (2010). 
[21] T. Brixner, J. Stenger, H. M. Vaswani, M. Cho, R. E. Blankenship, and G. R. Fleming, Nature 434, 625 (2005). 
[22] S. Savikhin, D. Buck, W. Struve, Chemical physics, 223, 303, (1997). 
[23] R. J. Sension, Nature 446, 740 (2007). 

[24] N. Renaud, M. A. Ratner and V. Mujica, J. Chem. Phys. 135, 075102 (2011). 
[25] S. M. Vlaming and R. J. Silbey, J. Chem. Phys. 136, 055102 (2012). 
[26] J. J. Hopfield, Proc. Nat. Acad. Sci. U.S.A., 71, 3640 (1974). 
[27] Th. Forster, Ann. Phys. 437, 55 (1948). 

[28] M. Mohseni, P. Rebentrost, S. Lloyd and A. Aspuru-Guzik, J. Chem. Phys. 129, 174106 (2008). 

[29] A. Olaya-Castro, Chiu Fan Lee, F. F. Olsen, and N. F. Johnson, Phys. Rev. B 78, 085115 (2008). 

[30] E. N. Zimanyi and R. J. Silbey, J. Chem. Phys. 133, 144107 (2010). 

[31] P. K. Ghosh, A. Y. Smirnov and F. Nori, J. Chem. Phys. 134, 244103 (2011). 

[32] M. Sarovar, A. Ishizaki, G. R. Fleming and K. B. Whaley, K. B., Nature Physics 6, 462 (2010). 

[33] A. Ishizaki and G. R. Fleming, J. Chem. Phys. 130, 234111 (2009). 

[34] D. Abramavicius and S. Mukamel, J. Chem. Phys. 133, 064510 ( 2010). 

[35] O. Miilken and T. Schmid, Phys. Rev. E 82, 042104 (2010). 

[36] M. B. Plenio and S. F. Huelga, New J. Phys. 10, 113019 (2008). 

[37] A. Ishizaki and G. R. Fleming, Proc. Natl. Acad. Sci. USA 106, 107255 (2009). 

[38] P. Rebentrost, M. Mohseni, I. Kassal, S. Lloyd, and A. Aspuru-Guzik, New J. Phys. 11, 033003 (2009). 
[39] F. Caruso, A. W. Chin, A. Datta, S. F. Huelga, and M. B. Plenio, J. Chem. Phys. 131, 105106 (2009). 
[40] F. Caruso, A. W. Chin, A. Datta, S. F. Huelga, and M. B. Plenio, Phys. Rev. A 81, 062346 (2010). 
[41] J. Striimpfer and K. Schulten, J. Chem. Phys. 134, 095102 (2011). 
[42] X. Chen and R. J. Silbey, J. Phys. Chem. B, 132, 205101 (2010). 

[43] M. A. Nielsen and I. L. Chuang, Quantum computation and information (Cambridge University Press, Cambridge, UK, 
2000). 

[44] R. Horodecki, P. Horodecki, M. Horodecki, and K. Horodecki, Rev. Mod. Phys.81, 865 (2009). 
[45] M. Horodecki, P. Horodecki and R. Horodecki, Phys. Lett. A 223, 1 (1996). 

[46] C. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. Peres, and W. Wootters, Phys. Rev. Lett. 70, 1895 (1993). 
[47] M. Horodecki, P. Horodecki, and R. Horodecki, Phys. Rev. A 60, 1888 (1999). 
[48] S. Popescu, Phys. Rev. Lett. 72, 797 (1994). 

[49] C.H. Bennett and S.J. Wiesner, Phys. Rev. Lett. 69, 2881 (1992) . 

[50] M. Hillery, V. Buzek, and A. Berthiaume, Phys. Rev. A 59, 1829 (1999). 

[51] C.A. Fuchs, N. Gisin, R.B. Griffiths, C-S. Niu, and A. Peres, Phys. Rev. A 56 1163 (1997) . 

[52] C. H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, A. Peres, W. K. Wootters, Phys. Rev. Lett. 70 1895 (1993). 

[53] V. Coffman, J. Kundu and W. K. Wootters, Phys. Rev. A 61, 052306 (2000). 

[54] H. Ollivier and W. H. Zurek, Phys. Rev. Lett. 88, 017901 (2001). 

[55] L. Henderson and V. Vedral, J. Phys. A 34, 6899 (2001). 

[56] V. Vedral, Phys. Rev. Lett. 90, 050401 (2003). 

[57] A. P. Majtey, A.R. Plastino and A. Plastino, "New features of quantum discord uncovered by q-entropies" arXiv:1112.4957 
[58] A. S. Davydov, Theory of Molecular Excitons (Plenum, New York, 1971). 

[59] D. P. Craig and S. H. Walmsley, Excitons in Molecular Crystals (Benjamin Inc., New York, 1968). 

[60] Y. Toyozawa, Optical Processes in Solids (Cambridge, New York, 2003). 

[61] A. Suna, Phys. Rev. 135, Alll (1964). 

[62] A. Thilagam and J. Singh, Phys. Rev. B 48, 4636 (1993). 

[63] A. Thilagam, Phys. Rev. A 81, 032309 (2010). 



16 



[64] T. Yu and J. H. Eberly, Phys. Rev. Lett. 93, 140404 (2004). 

[65] O. Giihne, F. Bodoky, and M. Blaauboer Phys. Rev. A 78, 060301(R) (2008). 

[66] A. Borras, A. P. Majtey, A. R. Plastino, M. Casas, and A. Plastino, Phys. Rev. A 79, 022108 (2009). 
[67] W. Diir and H. J. Briegel, Phys. Rev. Lett. 92, 180403 (2004). 

[68] L. Aolita, R. Chaves, D. Cavalcanti, A. Acm, and L. Davidovich, Phys. Rev. Lett. 100, 080501 (2008). 
[69] Zhong-Xiao Man, Yun-Jie Xia, and Nguyen Ba An, Phys. Rev. A 78, 064301 (2008). 
[70] O Giihne , F. Bodoky and M. Blaauboer, Phys. Rev. A 78 R060301 (2008). 

[71] A. Borras, A. P. Majtey, A. R. Plastino, M. Casas and A. Plastino, Phys. Rev. A 79 022108 (2009). 
[72] M. Hotta, Phys. Rev. A 80, 042323 (2009). 

[73] E. D Hondt and P. Panangaden, Quantum Inf. Comput. 6, 173 (2006). 
[74] M. Hillery, V. Buzek, and A. Berthiaume, Phys.Rev.A 59, 1829 (1999). 

[75] D. M. Greenberger, M. A. Home, A. Shimony, and A. Zeilinger, Am. J. Phys. 58, 1131 (1990). 
[76] A. R. R. Carvalho, F. Mintert, and A. Buchleitner, Phys. Rev. Lett. 93, 230501 (2004). 
[77] C. Simon and J. Kempe, Phys. Rev. A 65, 052327 (2002). 
[78] R. H. Dicke, Phys. Rev. 93, 99 (1954). 

[79] W. Diir, G. Vidal, and J. I. Cirac, Phys. Rev. A 62, 062314 (2000). 
[80] A. Kitaev and J. Preskill, Phys. Rev. Lett. 96 110404 (2006). 

[81] J. Batle, M. Casas, A. R. Plastino and A. Plastino, Phys. Lett. A, 296 251 (2002). 

[82] G. Vidal, R.F. Werner, Phys. Rev. A 65, 032314 (2002). 

[83] D. A. Meyer and N. R. Wallach, J. Math. Phys. 43 4273 (2002). 

[84] R. Chaves and L. Davidovich, Phys. Rev. A 82, 052308 (2010). 

[85] A. P. Majtey, A. Borras, A. R. Plastino, M. Casas and A. Plastino, Int. J. of Quantum Information 8, 505 (2010). 
[86] S. I. E. Vulto, M. A. de Baat, S. Neerken, F. R. Nowak, H. van Amerongen, J. Amesz, and T. J. Aartsma, J. Phys. Chem. 
B 103, 8153 (1999). 

[87] R. J. W. Louwe, J. Vrieze, A. J. Hoff, T . J. Aartsma, J. Phys. Chem. 101, 11280 (1997). 

[88] S. I. E Vulto, M. A. de Baat, R. J. W. Louwe, H. P. Permentier, T. Neef, M. Miller, H. van Amerongen, T. J. Aartsma, 

J. Phys. Chem. 102, 9577 (1998). 
[89] A. Thilagam, J. Chem. Phys. 136, 065104 (2012). 
[90] A. Thilagam, J. Phys. A: Math. Theor. 43, 155301 (2010). 

[91] H. P. Breuer and F. Petruccione, The theory of Open Quantum systems (Oxford, Oxford University Press, 2002). 
[92] Z. He, J. Zou, B. Shao and S. Y. Kong, J. Phys. B: At. Mol. Opt. Phys. 43, 115503 (2010). 
[93] A. Thilagam, J. Phys. A: Math. Theor. 42, 335301 (2009). 

[94] M. K. Sener, J D Olsen, C N Hunter, K. Schulten, Proc. Natl. Acad. Sci. USA 104, 15723 (2007). 
[95] D. Segale and V. A. Apkarian, J. Chem. Phys. 135, 024203 (2011). 



